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ABSTRACT 

We present results from photoionization models of low-metallicity Hn regions. These 
nebulae form the basis for measuring the primordial helium abundance. Our models 
show that the helium ionization correction factor (ICF) can be non-negligible for nebulae 
excited by stars with effective temperatures larger than 40,000 K. Furthermore, we find 
that when the effective temperature rises to above 45,000 K, the ICF can be significantly 
negative. This result is independent of the choice of stellar atmosphere. However, if 
an Hn region has an [O III] A5007/[O I] A6300 ratio greater than 300, then our models 
show that, regardless of its metallicity, it will have a negligibly small ICF. A similar, 
but metallicity dependent, result was found using the [O III] A5007/H/3 ratio. These 
two results can be used as selection criteria to remove nebulae with potentially non- 
negligible ICFs. Using our metallicity independent criterion on the data of Izotov & 
Thuan (1998) results in a 20% reduction of the rms scatter about the best fit Y — Z 
line. A fit to the selected data results in a slight increase of the value of the primordial 
helium abundance. 

Subject headings: galaxies: abundances — galaxies: ISM — Hn regions — ISM: abun- 
dances 



1. INTRODUCTION 

An accurate measurement of the primordial helium abundance would be an important test of 
standard big bang nucleosynthesis (Olive, Steigman & Skillman 1997), and would also constrain 
the values of the photon-to-baryon ratio and (Olive, Steigman & Walker 1999). The traditional 
procedure to measure the primordial helium abundance is to make use of the correlation between 
the helium mass fraction (Y) and metal abundance (Z). This correlation is then extrapolated to 
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zero metallicity to estimate the primordial mass fraction of helium, Y p . Spectroscopic observations 
of bright, low-metallicity extragalactic Hn regions provide the data for these studies (e.g., Olive & 
Steigman 1995; Olive, Steigman & Skillman 1997; Izotov, Thuan & Lipovetsky 1994, 1997; Izotov 
& Thuan 1998; Torres-Peimbert, Peimbert & Fierro 1989; Skillman, Terlevich & Terlevich 1998). 

To be cosmologically useful the value of Y p has to be determined to better than 5%. For- 
tunately, abundance determinations from measurements of line ratios is relatively straightforward 
(Peimbert 1975; Benjamin, Skillman & Smits 1999), and can, in theory, give the desired accuracy. 
However, to reach the needed level of precision, any systematic errors involved with target selec- 
tion, observations, and data analysis must be identified and corrected. Many such systematic errors 
have already been identified (Davidson & Kinman 1985; Dinerstein & Shields 1986; Pagel et al. 
1992; Skillman et al. 1994; Peimbert 1996; Izotov, Thuan & Lipovetsky 1997; Steigman, Viegas & 
Gruenwald 1997; Skillman, Terlevich & Terlevich 1998), but any errors resulting from the so-called 
"ionization correction factor" (ICF) have so far been assumed to be small. The ICF corrects for 
the fact that some amount of atomic (i.e., unseen) helium might be present in ionized regions of 
hydrogen (Osterbrock 1989; Peimbert 1975). This correction traditionally has been assumed to 
be zero because measurements of the primordial helium abundance employ observations of bright 
extragalactic Hn regions. These regions are excited by clusters of young stars with effective temper- 
atures greater than 40,000 K. Calculations by Stasinska (1990) and Pagel et al. (1992) showed that 
the helium ICF should be negligibly small for these Hn regions. As a result, recent determinations 
of Y p have assumed that the helium ICF is small. 

Very recently, Armour et al. (1999) presented calculations that showed that Hn regions excited 
by stars with temperatures greater than 40,000 K can have non-negligible ICFs. Armour et al. 
(1999) found that the ICFs were often negative (i.e., the helium ionized zone is larger than the 
hydrogen one; Stasinska (1980, 1982), Peha (1986)) for the hardest stellar continua. These results 
were confirmed by Viegas, Gruenwald & Steigman (2000). In this paper, we follow up on the work 
of Armour et al. (1999), and develop observational diagnostics of when the He ICF is important 
and when it can be ignored. We then apply these diagnostics to the data of Izotov & Thuan (1998) 
to illustrate how our technique can improve the precision of the measurement of Y p . 

We describe our calculations in § 2, and our results in § 3. The main results are summarized 
in S 4. 



2. DESCRIPTION OF CALCULATIONS 

In order to investigate the effects of a non-negligible ICF on the determination of the primordial 
helium abundance, we ran photoionization models of Hn regions and extracted the ICF for each 
nebula. These calculations are very similar to ones presented by Armour et al. (1999) and Bottorff 
et al. (1998), and were made with the development version of Cloudy, last described by Ferland et 
al. (1998). 
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Since we are modeling Hn regions, our models use the ISM abundances and grain model that 
were used and described by Armour et al. (1999). However, we scaled both the metal and grain 
abundances to lower values because, in this case, we are most interested in lower metallicity nebulae. 
The scaling was implemented so that all metals and grains were varied together relative to hydrogen 
and helium, but the He/H ratio was held constant. We modeled nebulae at three different metal- 
licities: 0/H=32, 64, and 128 parts per million [ppm] [Z = Z Q /23, Z /12, and Z Q /6 respectively). 
For each metallicity, 1936 models were computed for each of the following spectral energy distribu- 
tions: the LTE plane-parallel atmospheres of Kurucz (1991), the non-LTE, wind-blanketed, solar 
abundance CoStar atmospheres of Schaerer et al. (1996a, b), the earlier NLTE atmospheres of Mi- 
halas (1972), and for completeness blackbodies. We also ran models using the subsolar-abundance 
CoStar atmospheres which are spectrally slightly harder than the solar abundance ones. These 
models resulted in slightly more negative ICFs, but the values of the line ratio cutoffs (§ 3.2 k, 
§ 3.3) were not changed from the ones calculated with a solar abundance atmosphere. 

For each atmosphere we computed models with 10 cm -3 < n# < 10 6 cm" - 3 , 1(T 4 < U < 
KT - 25 , and 40,000 K< T eff < 50,000 K, where U is the ionization parameter defined as in Eq. 4 
of Armour et al. (1999). Giant extragalactic Hn regions that are observed are generally excited by 
large clusters, so this range of parameters should cover all such nebulae. We modeled the nebulae as 
plane-parallel constant density slabs, a simple way to characterize blister Hn regions. Our proposed 
diagnostic indicators are the [O III] A5007 and [O I] A6300 lines (§ 3.2 & § 3.3), and these should 
be fairly independent of the assumed geometry (sphere, sheet, or evaporating blister). The [O 
III] A5007/H/3 ratio represents the cooling per recombination (the Stoy ratio) and so is primarily 
sensitive to the stellar temperature (Stoy 1933; Kaler 1978) rather than geometry. Similarly, the [O 
I] A6300 /H/3 intensity ratio mostly measures the " softness" of the hydrogen ionization front (where 
the line forms - Netzer & Davidson 1979). 



3. RESULTS 

3.1. Temperature Dependence 

Figure 1 plots the ICF calculated from the Kurucz and CoStar models versus the stellar 
temperature. Note that we define the ICF such that an ICF of zero corresponds to zero correction: 

TCr ( H+ / H ) i m 
ICF " <He+/He) " (1) 

where the angle brackets denote the volume mean ionization fraction. This definition takes into 
account the presence of any He +2 in the nebula. Figure 1 clearly shows that one can obtain a non- 
negligible ICF for stars with temperatures greater than 40,000 K. The harder CoStar atmospheres 
give preferentially negative ICFs, which, if not taken into account, would result in a overestimate 
of the helium abundance. The same is true of the Mihalas atmospheres (not shown). These results 
agree with the calculations of Armour et al. (1999) and Viegas, Gruenwald & Steigman (2000). 
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The softer Kurucz atmospheres result in preferentially positive ICFs (i.e., the helium ionized zone 
is smaller than the hydrogen one), although, at temperatures greater than 45,000 K, they can also 
give negative ICFs. The blackbody atmospheres, the least realistic, are softer still, but even they 
can produce negative ICFs in some models at the highest temperatures. Therefore, negative ICFs 
seem to be found at high stellar temperatures independent of the type of stellar atmosphere. 

3.2. Metallicity Dependent Cutoff Criterion 

A negative ICF occurs as the results of penetrating high-energy photons preferentially ioniz- 
ing helium, due to its large photoionization cross section. This tends to be important for lower 
ionization parameter models, since these have significant regions where H and He are partially ion- 
ized. We expect these nebula to be characterized by lower [O III] A5007/H/? ratios (a measure of 
excitation) and larger [O I] A6300/H/3 ratios (since [O I] A6300 is formed in warm atomic regions). 

The results presented in § 3.1 show that it is not appropriate to simply assume that the ICF 
is zero when a nebula is excited by a star with a temperature greater than 40,000 K. However, it 
would be important to develop observational diagnostics for when the ICF is important and when 
it is not. Figure 2 shows such a diagnostic. A plot of ICF versus [O III] A5007/H/3 shows that 
beyond a line ratio of about 3-4 the ICF is negligible (for clarity results are shown for Kurucz and 
CoStar atmospheres only; a plot for blackbody and Mihalas atmospheres is very similar). However, 
the value of the cutoff will depend on metallicity. We found very small ICFs for line ratios greater 
than the following cutoff: 

([O III] A5007/H/?) Cutoff = (0.025 ± 0.004)O/H + (1.139 ± 0.306), (2) 

where O/H is measured in parts per million. Hn regions which have an [O III] A5007/H/3 ratio less 
than the cutoff for their metallicity might be subject to an ICF correction. Unless this correction 
can be made (and, in general, it cannot) these Hn regions should be removed from the abundance 
analysis, as they will increase the scatter in the dY/dZ relation used to determine the primordial 
helium abundance. 



3.3. Metallicity Independent Cutoff Criterion 

One can improve the above result by finding emission line ratios that should be independent 
of metallicity. Figure 3 plots the helium ICF versus the [O III] A5007/[O I] A6300 ratio. The figure 
shows that the ICF is negligible for a line ratio greater than about 300. Not surprisingly, this result 
is independent of metallicity. 

There are some Kurucz and blackbody models which result in a non-negligible ICF at line 
ratios larger than this cutoff. These models had low stellar temperatures (40,000-42,000 K), and 
were found over a narrow range in both \ogU (—1.5 to —2.25) and log njj (1.0-3.0). Their positive 
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ICF is a result of combining the low stellar temperatures and the softness of their atmospheres; their 
large [O III] A5007/[O I] A6300 ratio is a result of combining the fairly high ionization parameter 
with the low density. CoStar and Mihalas models, which are considered more "realistic", result in 
only very small ICFs in this region. 

Therefore, we find that any Hn region, regardless of its metallicity, that has an [O III] A5007/[O 
I] A6300 ratio less than about 300 might be subject to an ICF correction and should not be used 
to determine the primordial helium abundance. 

3.4. Application to Real Data 

To see how these new results affect the determination of Y p , we applied the metallicity in- 
dependent rejection criteria to the data of Izotov & Thuan (1998), and the results are shown in 
Figure 4. There are a number of points to note from this Figure: 

1. Our criterion rejects points over the entire range of metallicity, so there is no metallicity bias. 
The rejected points also fall evenly over the range of Y values, which implies that there is no 
correlation between Y and the [O III] A5007/[O I] A6300 ratio. 

2. At an [O III] A5007/[O I] A6300 cutoff of 300 there is a 20% reduction in the weighted rms 
scatter about the best fit line. This is encouraging evidence that part of the scatter was due 
to the ICF, and the situation has indeed improved by the implementation of the cutoff. 

3. The negative slopes predicted by using a cutoff > 300 result from weighted fits to the small 
number of data that remain after applying the cutoff, and are probably not realistic. Given 
the current data, a cutoff value greater than 300 is too severe (not practical). 

4. Despite the increased errors in the slope and the intercept, we find that implementing our 
cutoff will result in a larger value of Y p . Specifically, we find Y p = 0.2489 ± 0.0030 when 
the cutoff is taken at a [O III] A5007/[O I] A6300 ratio of 300. Izotov & Thuan (1998) find 
0.2443 ± 0.0015 (our fit gives 0.2443 ± 0.0013 with no cutoff), so the two results are barely 
consistent at the la level. This result is in the opposite direction of the shift predicted by 
the Monte Carlo simulations of Viegas, Gruenwald & Steigman (2000), but moves the value 
of Y p closer to the theoretically predicted values (Olive & Steigman 1995). 

There is always the possibility that systematic errors are introduced whenever data are rejected 
by a certain criterion. The above selection criterion preferentially selects Hn regions that have large 
[O III] A5007/[O I] A6300 ratios. This ratio is generally large whenever the [O III] A5007 line is 
strong (this is consistent with our metallicity dependent criterion). The [O III] A5007 line is a major 
source of cooling in a nebula, and so by selecting Hn regions with strong [O III] A5007 lines we are 
selecting regions ionized by hotter stars. Since these large extragalactic Hn regions are generally 
ionized by clusters, ones with strong [O III] A5007 lines are preferentially younger. However, this 
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is unlikely to introduce a systematic error in the primordial helium abundance determination, as 
young clusters can form at any metallicity. Indeed, Figure 4 shows that the points rejected by 
applying the cutoff span the entire range of metallicity. 

There is also the possibility that physical conditions within the Hn regions may bias our results. 
For example, in some nebulae the intensity of the [O I] A6300 line could be enhanced due to shock 
heating. This would lower the measured [O III] A5007/[O I] A6300 ratio, and could move it below 
our selection criterion. However, shock heating would not change the ICF of the nebula, so even 
though application of our criterion might reject such a Hn region, it will not bias the determination 
ofY p . 

Another potential situation in Hn regions is that the nebula may be matter bounded (i.e., 
optically thin to the Lyman continuum) in certain solid angles or sectors. Because there will be no 
hydrogen ionization front in such sectors, and therefore no [O I] A6300 line emission, the presence 
of such sectors will increase the measured [O III] A5007/[O I] A6300 ratio, possibly pushing it above 
our selection criterion. However, these sectors will also have little or no neutral helium within 
them, and so will have no ICF. Therefore, although the other, ionization bounded, sectors of the 
Hn region could have a non-negligible ICF, this will be diluted by the matter bounded sectors. We 
anticipate that even if such an Hn region were shifted into our selected data, there ought not be a 
large effect on determining Y p . According to Fig. 3, to severely bias the results a number of points 
would have to shift rightward in [O III] A5007/[O I] A6300 by a factor larger than ten; but such 
a large shift in the line ratio would probably result in a large dilution in the ICF. More modeling 
would be needed to quantify how little an impact matter bounded sectors would have on the [O 
III] A5007/[O I] A6300 ratio and ICFs. 



4. CONCLUSIONS 

In this paper we have shown the following: 

1. There can be a non- negligible ICF correction for Hn regions excited by stars with tempera- 
tures greater than 40,000 K. At temperatures higher than 45,000 K, the ICF is preferentially 
negative. This result is independent of the atmosphere of the O star. 

2. There is a simple procedure to determine if an ICF correction needs to be made for a given Hn 
region. If the [O III] A5007/[O I] A6300 ratio is greater than 300, then no correction is needed. 
This criterion is independent of metallicity. If the [O I] A6300 line cannot be measured, then 
there is a metallicity dependent cutoff (Eq. 2) that can be used with the [O III] A5007 line. 

3. Applying the metallicity independent criterion to the data of Izotov & Thuan (1998) results 
in reducing the rms scatter about the best fit Y — Z line by 20%. This will help remove 
systematic errors relating to unrecognized ICF effects, and ought to improve the reliability 
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of the Y p determination. Furthermore, an analysis of the selected data gives a larger value of 
Y p than was originally measured, which is closer to the theoretically expected value. 
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Fig. 1. — A plot of the He ICF values obtained from the CoStar and Kurucz photoionization 
models versus stellar temperature of the atmospheres. Only points with ICF between ±10% are 
plotted. There are non-negligible ICFs at all values of stellar temperature. The CoStar atmo- 
spheres generally result in negative ICFs, while the softer Kurucz atmospheres give positive ICFs 
for temperatures less than about 45,000 K. These models were run with a metallicity of (0/H)=64 
ppm. 

Fig. 2. — A plot of He ICF vs. [O III] A5007/H/3 using data from the CoStar and Kurucz grids 
run at (0/H)=64 ppm. Only points with ICF between ±10% were plotted. At this metallicity the 
ICF-related cutoff is at a relative line strength of about 3-4. See the text for discussion on how 
cutoff varies with metallicity. 

Fig. 3.— a) A plot of He ICF vs. [O III] A5007/[O I] A6300 with data obtained from the CoStar 
and Kurucz model grids run at (0/H)=64 ppm. The ideal ICF-related cutoff here is at a line ratio 
of 300. This cutoff is independent of metallicity. The Kurucz models that have positive ICFs in the 
allowed zone are models with a particular combination of parameters (see text), b) Like a) but for 
Mihalas and blackbody atmospheres. 

Fig. 4. — This figure shows the effects of applying the metallicity independent rejection criterion 
to the data of Izotov & Thuan (1998). The solid line is a weighted least-squares fit to the selected 
data points shown by the solid symbols. The open symbols are the rejected points. For reference, 
the dashed line is the fit with no cutoff applied. Note that as the [O III] A5007/[O I] A6300 cutoff 
becomes larger, the scatter of the points about the best-fit line becomes smaller (there is a 20% 
reduction in the rms scatter when the cutoff is 300). This procedure also shows that the value of 
Y p determined by Izotov Sz Thuan (1998) might be an underestimate. 
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